Entanglement generation outside a Schwarzschild black hole and 

the Hawking effect 

Jiawei Hu 1 and Hongwei Yu 1 ' 2 * 
1 Institute of Physics and Key Laboratory of Low Dimensional 
Quantum Structures and Quantum Control of Ministry of Education, 
Hunan Normal University, Changsha, Hunan £10081, China 
2 Center for Nonlinear Science and Department of Physics, 
Ningbo University, Ningbo, Zhejiang 315211, China 

Abstract 

We examine the Hawking effect by studying the asymptotic entanglement of two mutually in- 
dependent two-level atoms placed at a fixed radial distance outside a Schwarzschild black hole in 
the framework of open quantum systems. We treat the two-atom system as an open quantum 
system in a bath of fluctuating quantized massless scalar fields in vacuum and calculate the con- 
currence, a measurement of entanglement, of the equilibrium state of the system at large times, 
for the Unruh, Hartle-Hawking and Boulware vacua respectively. We find, for all three vacuum 
cases, that the atoms turn out to be entangled even if they are initially in a separable state as long 
as the system is not placed right at the even horizon. Remarkably, only in the Unruh vacuum, 
will the asymptotic entanglement be affected by the backscattering of the thermal radiation off the 
space-time curvature. The effect of the back scatterings on the asymptotic entanglement cancels 
in the Hartle-Hawking vacuum case. 
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I. INTRODUCTION 



Classically, black holes are described as massive objects with such a strong gravitational 
field that even light cannot escape from them. However, Hawking finds, in the framework 
of quantum field theory in curved space-time, that a black hole is not completely black, but 
emits thermal radiation with a black body spectrum due to quantum effects Ever since 
this surprising discovery, the Hawking effect has attracted widespread interest in the physics 
community and extensive works have been done trying to understand it in various different 
physical contexts (See for example, Ref. [irul). 

In this paper, we will try to examine the Hawking effect in terms of the entanglement 
generation in the framework of open quantum systems. The open system we are going to 
study consists of two mutually independent static two-level atoms subjected to a bath of 
fluctuating quantized massless scalar fields in vacuum outside a Schwarzschild black hole. 
We will analyze the time evolution of the density matrix describing the system using the 
well-known techniques in the theory of open quantum systems. Let us note that the reduced 
dynamics of a single static detector (a two-level atom) placed outside a Schwarzschild black 
hole interacting with quantized massless scalar fields in the Unruh, Hartle-Hawking and 
Boulware vacua has been recently investigated [9( and it has been found that in both the 
Unruh and Hartle-Hawking vacua, the detector will spontaneously excite with a nonvanishing 
probability as if there were thermal radiation at the Hawking temperature. An another way 
to look at the issue is to study the equilibrium state of the detector. In this regard, it has 
been shown that the detector is asymptotically driven to a thermal state at the Hawking 
temperature in the spatial asymptotic region, regardless of its initial state [lO]. This approach 
has also been applied to reproduce both the Unruh effect and the Gibbons- Hawking effect, 
in Ref. [U[ and [l2j, respectively. 

If the single atom is replaced by two independent two-level atoms, the situation becomes 
physically more interesting. A study of the reduced dynamics of the two-atom system may 
reveal whether the asymptotic equilibrium state of the system will be entangled or not. It 
is known that an environment usually leads to decoherence and noise, which may cause 
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entanglement that might have been created before to disappear. However, in certain cir- 

-lisj. 



cumstances, the environment may enhance entanglement rather than destroying it |13 
The reason is that an external environment can also provide an indirect interaction between 
otherwise totally uncoupled subsystems through correlations that exist. It is of interest to 
see whether a bath of fluctuating vacuum scalar fields outside a Schwarzschild black hole 
can provide such an indirect coupling to enhance the entanglement. This is what we are 
going to pursue in the present paper. We will, in the hope of gaining new understanding of 
the Hawking effect from a different perspective, examine the asymptotic entanglement of the 



two-atom system in a bath of massless scalar fields in the Unruh [19] . Hartle-Hawking [20], 
and Boulware 21] vacua outside a black hole. Let us note here that a similar issue related 



to the Unruh effect has already been worked out by Benatti and Floreanini, in which they 
have studied a uniformly accelerating two-atom system 11] and found that the two initially 



separable atoms will become entangled, just as if they were immersed in a thermal bath at 
the Unruh temperature. 



II. THE MASTER EQUATION 

The system we study consists of two mutually independent static two-level atoms of van- 
ishing separation in interaction with a bath of fluctuating quantum scalar fields in vacuum 
outside a Schwarzschild black hole which is described by the following metric 

ds 2 ={ l - 2J f) - Y^Mfr ~ ^ + Sin2 M<p2) ■ (1) 

Although the atoms are mutually independent, the fluctuating vacuum fields with which the 
atoms are coupled may provide an indirect interaction and therefore a means to generate 
entanglement between them. In this paper, we are particularly concerned with the issue 
of whether the atoms can be entangled when system reaches an equilibrium state as time 
becomes large. Generically, the total Hamiltonian of the system takes the form 

H = H S + H^ + XH' . (2) 
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Here H* is the Hamiltonian of the two atoms, 



H s = hU+hW, H^ = ^rn<r\ a \ (a = 1,2), (3) 

where = oi <g> oq, <x- 2 ' ) = a <g> cr^, <Tj (z = 1, 2, 3) are the Pauli matrices, <7o the 2x2 unit 
matrix, n = (ni, n 2 , n 3 ) a unit vector, and Uq the energy level spacing. is the Hamiltonian 
of free massless scalar fields. The Hamiltonian H' that describes the interaction between the 
two atoms with the external scalar fields is given by 

3 

H ' = J2 ® a o)^(t, xi) + (a (8) a^^t, x 2 ) ] . (4) 
We assume that the scalar fields can be expanded as 

TV 

= E K4> { -\x) + (x»r4> {+ \x)} , (5) 

a=l 

where <j>^{x) are positive and negative energy field operators of the massless scalar field, 
and Xu are t ne corresponding complex coefficients. 

At the beginning, the whole system is characterized by the total density matrix p to t = 
p(0) g) |0)(0|, in which p(0) is the initial reduced density matrix of the two-atom system, and 
|0) is the vacuum state of field $(x). In the frame of the two-atom system, the evolution in 
the proper time r of the total density matrix p tot satisfies 

= -iL„Mr)] , (6) 

where the symbol Lh represents the Liouville operator associated with H 

L H [S] = [H,S}. (7) 

We assume that the interaction between the atoms and the field is weak, i.e., the coupling 
constant A in (J2]) is small. In the limit of weak coupling, the evolution of the reduced density 
matrix p(r) can be written in the Kossakowski-Lindblad form [ll, 17, 22], 23] 

= -i[H eS , p(r)] + C[p(r)) , (8) 
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with 

H« = H.- t -j2Y,B ij a!r ) ay>, (9) 

a, 0=1 i,j=l 

and 

£ W ^tt C^afpaf - afaf „ - pafaf] . (10) 

a,P=l i,j=l 

The coefficients of the matrix Cy and are determined by the Fourier and Hilbert trans- 
forms of the field correlation functions 

Gy(x-y) = <0|* 4 (a;)* i (y)|0> , (11) 

which are defined as follows 

/>oo 

G l] {X)= / dre iXT G l3 (T) , (12) 
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oo " " J — oo 



/Cy(A) = / rfr sign(r) e** Gy (r) = ^ [°° du ^4 , ( 1 3 ) 



a; — A 



in which P denotes principal value. It can be shown that the Kossakowski matrix Cy can 
be written explicitly as 

c *i= E E^(^)4?W;°, (14) 

£=+ ,-,0*,J=l 

where 

= Tii % , = i - Tij ± iey fc n fc ) . (15) 

Similarly, the coefficients of can be obtained by replacing Gki{&) with K.m(£uj) in the 
above equations. 

If we assume that the field components are independent, or the coefficients satisfy 

N 

E >W = V , (16) 

a=l 

the field correlation functions in fllip are diagonal such that 



Gij(x - y) = SijG(x - y) , (17) 
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in which G(x — y ) is the standard Wightman function, and the Kossakowski matrix Cy can 



be written as 



111 



Cij = A 5ij - iB e ijk n k + C UiUj , (18) 



where 



a = l -\g{u ) + g(-u )} , b = l -\g{ujv) - g{-u Q )] , c = g(o) - a . (19) 



III. THE ASYMPTOTIC ENTANGLEMENT 

In order to figure out whether the entanglement will be generated between the atoms and 
whether it can persist asymptotically through the analysis of the master equation governing 
the density matrix, we need to calculate the coefficients which are determined by the field 
correlation function in the vacuum state. However, when a vacuum state is concerned in a 
curved space-time, a delicate issue then arises as to how the vacuum state of the quantum 
fields is determined. Normally, a vacuum state is associated with non-occupation of positive 
frequency modes. However, the positive frequency of field modes is defined with respect 
to the time coordinate. Therefore, to define positive frequency, one has to first specify a 
definition of time. In a spherically symmetric black hole background, three different vacuum 



states, i.e., the Unruh [19], Hartle-Hawking [20|], and Boulware |2jJ vacuum states, have been 



defined, each corresponding to a different choice of time coordinate. In the next section, we 
examine the asymptotic entanglement for the atoms in all vacua. 



A. The Unruh vacuum 

Let us begin our discussion with the Unruh vacuum which is supposed to be the vacuum 
state best approximating the state following the gravitational collapse of a massive body 
to a black hole. In order to analyze the asymptotic state of the system, we will deal with 
the master equation ()8]) and see how the reduced density matrix of the two-atom system p 
evolves. 
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Generally, we can write the reduced density matrix of the two-atom system in the form of 
1 



P(r) , 



cr ® cr + Poi(T) ctq ® o"i + Pio(r)ai <g> a + Pij{r)oi <g> ct\,- 



(20) 



which is normalized as Tr(p) = 1 with det(p) > . Here we are interested in the final 
equilibrium density matrix p°°, which does not change with time, i.e., d T p°° = 0. Let us note 
that the unitary term involving Hamiltonian H e g in the master equation can be ignored since 
it does not give rise to the entanglement phenomena jllj] and we only need to examine the 
effects produced by the dissipative term C[p] in (JHD- As a result, the equilibrium condition 



becomes C[p°°] = 0. After some direct calculations, one obtains the components of p°° rjj 



where R = B/A, is the trace of the density matrix r* = S? =1 pjj(r), which is actually a 
constant of motion, and the positivity of p(0) requires — 3 < r* < 1. 

In order to determine whether the final equilibrium state is entangled or not, we take the 
concurrence as a measurement of the entanglement, which is defined as C[p] = max{0, Ai — 
A2 — A3 — A4} , where A M (p = 1, 2, 3, 4 ) are the square roots of the non-negative eigenvalues 
of the matrix p(a 2 <E> 02)p T (02 ®<J2) in decreasing order, and T stands for transposition. The 
value of C[p] ranges from 0, for separable states, to 1, for maximally entangled states. For 
the current case, the concurrence is 



C(p°°) = max 
which is non-zero provided obeys 



(3 - R 2 ) 
2{3 + R 2 ) 



5R 2 -3 
3-R 2 ~ T * 



, (22) 



5# 2 - 3 

r. < -j^ ■ (23) 



This result implies that as long as the condition ( 1231) is satisfied, the equilibrium state will 
turn out to be entangled, even if the initial state is separable. 

In order to make our discussion more concise, let us come to a simple example. The initial 
state of the system is taken to be a separable state provided by the direct product of two 



pure states: 



p(0) = p n ® p n 



p n = - (l + ft ■ <i) , p m — - ( 1 in ■ <r 



(24) 



2 V J ' 2 

where n and m are two unit vectors. Here r* = n • m, so the asymptotic entanglement is 
maximized when n — — in: 

C\P°°] = • (25) 



3 + R 2 

Eq. (|25p shows that, in our simple case, the concurrence increases with R monotonically 
when R ranges from to 1, and it reaches its maximum = 1/2 for R — 1. So, we need 

only evaluate R in order to analyze the entanglement of the system. 

Now let us begin to compute R in the Unruh vacuum case. The Wightman function for 
massless scalar fields in the Unruh vacuum is given by 24j-l26| 



G+(*,*') = E 



ml 



Y, 



lm\ 



i^,r)| 2 

I p — 2'KU) J K 



6{uj)\%(uj,r)\ 2 



(26) 



where k = 1/4M is the surface gravity of the black hole. Its Fourier transform is 



T 



1 oo p 

6{\^){l + 2l)\%{\^r)\ 2 

57TA * — ' 



-4 



+ 



{l + 2l)\R l {\ y /g^,r 



1=0 



where we have used the relation 



m=—l 



21 + 1 

47T 



(27) 



(2f 



and here K r is defined as ^/y^oo • The above Fourier transform, which is needed in our 
computation of R, is hard to evaluate, since we do not know the exact form of the radial 
functions Ri(u, r) and Ri(u>, r). Here, we choose to compute it both close to the event horizon 
and at infinity. For this purpose, let us recall that the radial functions have the following 
properties in asymptotic regions [261 ]: 

4u 2 



oo 

^(2Z + l)|5,(a;,r) 



I _ 2M ' 



oo 

-^(2/ + l)|^M 



r 2M 



r — > oo 



(29) 



1=0 



( 1 

ZM2 5> 



z=o 



«=0 



4M 2 
4w 2 



i _ 2M ' 
v 1 



1) I A (00 



r -> 2M 



r — t- oo 



Inserting Eq. ([27]) into Eq {H}, and using Eq. ([29]) and ([30]), one finds that 
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r -> 2M : 



(USD 


and 




W r 




4tt [ 




W r 




4tt [ 



and 



r — >■ oo : 



5 



^ [1 + (?oo /KVto, r) + y Poo/^oVW, 0] 

Wo 
47T 



[1+^00 /(woV^OO) r )] 



where f(u,r) is defined as 



/(w,r) 



1 



4r 2 cu 2 



£(2J + l)|B,(w) 



Straightforward calculations then yield in the asymptotic regions, 

' 1 + googa^, 2M) 

R _B _ 1 + 9oof(^oy/9oo, 2M) + e2TM0 / Kr _ 1 ' 
^4 I 1 + 5 , oo/(w 0% /^oo> r ) 



r — >■ oo 



(30) 



(31) 



(32) 



(33) 



(34) 



k 1 + 9oof(uo^goo, r) + ^77— y^oo/(w 0v /5oo, r) 

We can see that, in the vicinity of the event horizon, the first two terms in the denominator are 
the same as their counter-parts in the numerator, and the third one is the standard Planckian 
factor. At infinity, the Planckian factor is modified by a grey-body factor j7oo/(wo\/^oo> r ) 
caused by the backscattering off the space-time curvature. This suggests that at the horizon, 
there is a thermal flux going outwards, which is weakened by the backscattering off the 
curvature on its way to infinity. Actually, it has been shown that, in the framework of open 
quantum system, the spontaneous excitation rate per unit time of a particle detector from 
the initial ground state i to the final excited state / is just 



r^ f = 2(A -B) = 2g(-u ) 
9 



(35) 



This means that the difference between the denominator and the numerator of R is pro- 
portional to the spontaneous excitation rate per unit time of a particle detector, i.e., the 
strength of the thermal radiation. The larger the difference between A and B, the stronger 
the thermal radiation, and the less the two-atom system gets entangled. 

Here we note that, using the geometrical optics approximation 24j, the transmission 



amplitude Bi(uj) can be approximated as Bi(oo) ~ 6(y/27Muj — /), where 6(x) is the standard 
step function, which gives for x < and 1 for x > 0. So #oo/( w oa/9oo; r ) can ^ e simplified 

as 

27M 2 g 00 27M 2 ( 2M 



WfoVflto, r) « = _^l— -) S /(r) , (36) 

and in both the asymptotic regions, goo/( w oy(9oo> r ) — * 0. Allowing for this, we have, in the 
vicinity of the event horizon, 

e 2lTUl /K r _ 1 

R = - 9 — 7 . (37) 

It is interesting to note that the R we obtain here is exactly the same as that in the case of 
a two-atom system immersed in a thermal bath at the temperature T = K r /2ir = T# / y/goo 
in a flat space-time, where T# = n/2n is the Hawking temperature [1]. When r — )■ 2M, the 
temperature T is divergent, since this effective temperature is a result of both the thermal 
flux from the black hole and the Unruh effect due to that the system is accelerating with 
respect to the local free-falling inertial frame so as to maintain at a fixed distance from the 
black hole, and the acceleration diverges at the horizon. In this case, the concurrence 
approaches zero, which means that final equilibrium state of the two atom system, which 
is very close to the event horizon, will not be entangled. As the atoms are placed farther, 
then the thermal radiation becomes weaker due to the back scattering off the space-time 
curvature and the concurrence grows larger. At the infinity, the grey-body factor vanishes 
and R approaches 1, so that the concurrence C[p°°] tends to reach its maximum 1/2. This 
result suggests that, in the Unruh vacuum, no thermal radiation is felt at the infinity due to 
the back scattering of the outgoing thermal radiation off the space-time curvature. Here it 
is clear that in the vicinity of the horizon, the entanglement of the system is enhanced due 
to the back scattering, while in the infinity, the concurrence is smaller than its maximum if 
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we allow for the thermal radiation emitted from the horizon although it is back scattered by 
the space-time curvature. 

Compared with the case of a thermal bath in a flat space-time, we find that the parameters 
A and B are modified by the grey-body factor f(u, r), which is a function of M and r for a 
given energy gap ojq. A similar result is derived when studying the entanglement generation 
in atoms immersed in a thermal bath of external quantum scalar fields with a reflecting 
boundary j^]. This result implies that the back scattering of vacuum field modes off the 
space-time curvature of the black hole manifests in much the same way as the reflection of 
the field modes at the reflecting boundary in a flat space-time. 



B. The Hartle-Hawking vacuum 



Now let us move on to the Hartle-Hawking vacuum. The Wightman function is now given 



by 



24-26 



G + (s,s') = £ 



ml 



x l Y lm (e,4>)\ 2 



Amu 



jAt 



J g— 2ttoj/k 



Ri(u,r) z + 



%(tu,r)\ 2 



and its Fourier transform is 



duj , 



(3? 



e iXT G+{x,x')dr 
-4 



£ 

1=0 



;i+2Q 

8ttA 



^(Av/<7oo> r ) I 2 , I % 



+ 



~y g— 27tA/kv 



Similar calculations then lead to 



r^2M: 



A - ^ eWftr _ 1 [1 + Soo /Kv^, 2M)] 
B^^[l+g 00 f(u ^,2M)}, 



(39) 



(40) 



r — > oo : 



.4 
B 



u e 2 ™ o/Kr + 1 



wo 

47T 



- _ 1 [1 + #oo f(uoV9oo, r)] 



[l + ^oo /(w ov / <7oo>)] 



(41) 
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and 



l + ^oo/(wov^oO' 2M ) nnT 
- , r — > 2M , 



r — > oo . 



1 + g 00 f(u ^g^, 2M) + e2 ^ % r _ 1 + g^^ oo/C^oV^ 2M ) 
1 + 9oof(uoy/goo,r) 



1 + goof(uoy/goo, r) + e 2^ /K r _ 1 + e 2^ / KT . ^ goofiwo y/goo, r) 

(42) 

Unlike that in the Unruh vacuum case, here R is the same in both the asymptotic regions. 
We have two Planckian terms in the denominator. One is the standard one and the other 
is a Planckian factor modified by a grey-body factor caused by the backscattering off the 
space-time curvature. This suggests that there are thermal radiation outgoing from the 
horizon and that incoming from infinity, both of which are weakened by the backscattering 
off the curvature on their way. Therefore, the Hartle-Hawking vacuum is actually a state 
that describes a black hole in equilibrium with an infinite sea of blackbody radiation 26]. 
Here we notice that Eq. (I4"2j) can be simplified as 

R = ~, = ^T—i 7 ■ ( 43 ) 

At close to the horizon, i.e., when r — y 2M, R — y 0, which means that final state of the 
two-atom system is not entangled. When r — > oo, R — > ^/^| ) which is the maximal value 
it can reach in this case and it is nonzero, so two atoms will be entangled even if they are 
separable initially. Here, although both A and B are modified by a grey-body factor f(u,r), 
the modification cancels when we evaluate R. So, R is the same as what we get in the case 
of a thermal bath. This result shows that, both in the vicinity of the horizon and at infinity, 
the impact of the black hole to the equilibrium entanglement of the system is the same as 
that of a thermal bath at temperature T = T^/y^oo in a flat space-time. As r — > oo, the 
effective temperature becomes the Hawking temperature, since the acceleration needed to 
maintain the two-atom system at a fixed distance vanishes, and the temperature is purely 
due to the thermal bath the black hole immersed in. 
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C. The Boulware vacuum 



For the Boulware vacuum, the Wightman function is given by 24j-|26] 

/oo — iojAt , j 
1 Y lm (9, 0) | 2 [| Ri(uj, r) | 2 + | r) \ 2 ]du . (44) 

The Fourier transform with respect to the proper time is 

e^G + [x(r)]dr = V J?(Av5to,r) | 2 + | ^(A^r) | 2 ]^(A) , (45) 

mi 87rA 

where 0(A) is the step function. Plugging Eq. (145]) into Eq (|T9|) . we have 

A = 5 = El7i ±i n^(^0| 2 + lfe,r)| 2 ] . (46) 

So 



167TW 
«=0 



« = § = (47) 
everywhere outside the event horizon, and the concurrence of the equilibrium state is 

C[p°°] = C[p°°] max = \ (48) 

which is the same as that in a Minkowski vacuum. No thermal radiation is present. 



IV. CONCLUSION 

In summary, we have examined the asymptotic entanglement between two mutually inde- 
pendent two-level atoms at a fixed radial distance outside a Schwarzschild black hole in the 
paradigm of open quantum systems. We treat the two-atom system as an open quantum 
system in a bath of fluctuating quantized massless scalar fields in vacuum and have studied 
the Boulware, Unruh, and Hartle-Hawking vacua respectively. 

In the Hartle-Hawking and the Unruh vacuum cases, the concurrence attains a non-zero 
value less than 1/2 as long as the two-atom system is not placed right at the event horizon, 
indicating that the two atoms will turn out to be entangled even if they are initially separable. 
For the Unruh vacuum case, the concurrence is the same as if there were an outgoing thermal 
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flux of radiation from the event horizon, which is backscattered by the curvature of the 
space-time. For the Hartle-Hawking vacuum case, the concurrence behaves as if the atom 
were in a thermal bath of radiation at a proper temperature which reduces to the Hawking 
temperature in the spatial asymptotic region. Remarkably, only in the Unruh vacuum, will 
the asymptotic entanglement be affected by the backscattering of the thermal radiation off 
the space-time curvature. The effect of the back scatterings on the asymptotic entanglement 
cancels in the Hartle-Hawking vacuum case. 
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